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Integral  ProbabillW  of  Auroral  Electron  Flux  Events  from 

SSJ/4  DMSP  F9  Electron  Measurements 


1.  INTRODUCTION 


Determination  of  the  ^obal  pattern  of  auroral  particle  precipitation  is  important  both  for 
the  understanding  of  the  physical  processes  producing  and  maintaining  the  auroral  oval  and 
for  the  prediction  of  the  effects  of  auroral  particle  precipitation  on  operational  Air  Force 
systems.  For  this  reason,  an  extensive  set  of  studies  has  been  conducted  to  quantify  the  global 
pattern  using  measurements  from  the  SSJ/3  and  SSJ/4  particle  spectrometers  flown  on  the 
spacecraft  of  the  Defense  Meteorological  Satellite  Program  (DMSP)  1.2,3,4,5,6,7  These  studies 
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have  been  of  two  types:  1)  studies  nidiere  the  behavior  of  the  average  electron  and  ion  spectrum 
Is  determined  as  a  function  of  either  the  level  of  geomagnetic  activity  or  the  value  of  solar 
wind  parameters,  and  2)  studies  where  the  probability  of  occimence  of  different  levels  of 
auroral  precipitation  is  determined  as  a  function  of  the  level  of  geomagnetic  activity. 

The  DMSP  data  set  Is  well  suited  for  such  studies  because  of  Its  size.  The  SSJ/3  and  SSJ/4 
detectors  are  operated  continuously  on  the  DMSP  satellites,  recording  complete  electron  and 
Ion  spectra  once  per  second.  Approximate^  15  years  of  the  data  Is  available,  consisting  of 
approximate^  150  million  spectra.  Such  a  large  data  set  allows  determination  of  the  averages 
and  probability  distribution  to  a  high  precision. 

The  present  study  is  an  extension  of  this  work.  This  study  uses  the  1990  measurements 
from  the  F9  satellite  taken  In  the  1030-2230  meridian.  The  purpose  of  this  study  is  to 
determine  not  only  the  probability  of  observing  a  specified  level  of  auroral  precipitation  at  a 
given  spatial  location,  but  also  the  Integral  probability  of  observing  this  same  level  over  a 
range  in  latitude  starting  from  a  particular  location.  Only  the  first  part  of  these  Issues  was 
addressed  In  earlier  studies  of  DMSP  data  from  satellites  F2  and  P4.^  The  results  of  such  a 
study  allow  the  determination  of  the  Integral  probability  of  occurrence  for  auroral 
precipitation  events  above  any  specified  level  In  Integral  energy  flux  and  average  energy  for 
any  line  of  sight  through  the  auroral  region.  This  is  an  Important  tool  for  determining 
background  and  Interference  signals  for  any  system  operating  In  the  auroral  region. 


2.  ANALYSIS  TECHNIQUE 


The  purpose  of  this  study  Is  to  determine  the  probability  of  observing  auroral  electron 
precipitation  at  a  specified  Intensity,  si>ectral  hardness,  and  level  of  geomagnetic  activity  both 
at  a  specific  location  In  the  oval,  and  Integrated  over  ary  path  throu^  the  auroral  oval.  The 
data  for  this  study  are  from  the  SSJ/4  sensor  on  the  DMSP/F9  spacecraft.  All  data  taken  in 
1990  were  used  hi  the  study.  The  F9  satellite  Is  In  a  circular  sim  synchronous  oiblt  tn  the 
1030-2230  magnetic  meridian  at  an  altitude  of  approximately  840  kilometers.  The  spatial 
coverage  of  the  F9  orbit  In  a  corrected  geomagnetic  latitude  (CGL)  magnetic  local  time  (MLT) 
coordinate  system  is  shown  In  Figure  1 . 

For  this  stu^  we  constructed  a  series  of  spmtial  arrays.  In  each  such  array  we  divided  the 
high  latitude  region  Into  8  three-hour  sectors  in  MLT.  Each  array  covered  the  latitude  range 
fiom  50®  to  90®  CGL.  One  such  array  was  constructed  for  each  of  six  thresholds  In  integral 
energy  flux  and  average  energy  (spectral  hardness)  and  for  four  levels  of  geomagnetic  activity 
as  measured  by  Kp.  Thus  there  were  a  total  of  144  arrays  (6  thresholds  In  Integral  energy  flux 
times  6  thresholds  in  average  energy  times  4  I4>  ranges).  The  thresholds.  ranges,  and 
spatial  bins  used  are  listed  In  Table  1.  The  equivalent  radiance  bins  are  set  assuming  that  an 
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electron  energy  Input  of  1  erg/cm^-tiec  produces  1  Idlo-Raylel^  of  hindno^ty  at  a  wavdength 
ot  3969  A  [Figure  12. 19.  Chapter  12.  In  Reference  8 1. 


Table  1.  Input  Parameter  Ranges  and  Bins 


Parameter 

Bins 

Energy  flux  thresholds 

(1)  0.5  (2)  1.6  (3)  5.0  (4)  16.0  (5)  50.0  6}  160.0  (kR) 

(1)  108  (2)  108.5  (3)  io9  (4)  109-5  (5)  lOlO  (6)  1010  5 

KeV/(cn^sec-str) 

Average  energy 
thresholds  (KeV) 

(1)  1.0  (2)  3.0  (3)  7.0  (4)  10.0  (5)  13.0  (6)  17.0 

MLT  sector  (hours) 

(1)  0  -  3  (2)  3  -  6  (3)  6  -  9  (4)  9  -  12 

(5)12-15  (6)15-18  (7)18-21  (8)21-24 

I4>  ranges  or  midnight 
Equatorward 
Botmdaty 
(MEB  in  0) 

Kp  (1)  Oq  .  1+  (2)  2-  .3+  (3)  4-  .  5+  (4)  6-  .  9o 

Latitude  zones 

20  latitude  bins  from  50°  to  90°  with  2°  Increments 

In  each  MLT  sector  of  each  array  we  constructed  a  20  x  20  matrix.  The  diagonal  dements 
of  the  matrix  represent  the  probability  of  seeing  an  dectron  precipitation  event  above  the 
thresholds  for  that  array  during  a  pass  over  the  auroral  zone.  The  20  diagonal  dements  cover 
the  latitude  range  horn  50  to  90**  CGL  In  2**  steps.  For  a  given  row  In  the  matrix,  the  dements 
to  the  left  of  the  diagonal  represent  the  Integral  probability  of  seeing  an  event  above  the 
thresholds  for  that  array  Integrated  from  the  latitude  of  the  diagonal  to  a  latitude  equatorward 
of  that  latitude.  Each  Individual  position  to  the  left  represents  an  additional  Int^ratlon  over 
2°  of  latitude  equatorward.  Similarly,  the  dements  to  the  rl^t  represent  the  Integral 
probability  poleward  of  the  latitude  of  the  diagonal  element. 

Ihe  values  in  the  matrix  were  determined  as  follows;  For  each  auroral  pass,  the  Integral 
energy^  flux  and  average  energy  are  calculated  for  eadi  1-second  spectrum  returned  by  the 
SSJ/4  detector.  We  assume  that  the  pattern  of  auroral  dectron  precipitation  is  time 
stationary  over  the  period  of  the  auroral  pass.  Based  on  the  values  of  the  Integral  energy  flux 
and  average  energy  flux  calculated,  along  with  I4>  at  the  time  the  spectrum  was  taken,  we 
determined  to  which  If  any  of  the  arrays  the  event  was  to  be  assigned.  Based  upon  the  MLT. 
the  appropriate  matrix  was  chosen  and  for  that  matrix  the  CGL  Is  used  to  determine  the 
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Figure  1.  DMSP  F9  Satellite  Coverage  in  1990 
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ctHTCct  diagonal  element  to  adilch  the  event  applies.  If  the  spectrum  meastired  Is  the  first 
spectrum  for  that  pass  observed  for  that  set  of  thresholds  and  MLT  and  for  that  diagonal 
element,  then  we  Increment  the  counter  for  that  diagonal  element  by  one.  If  for  that  pass  a 
spectrum  for  that  latitude  and  MLT  has  already  been  ol»served  above  the  thresholds,  then  the 
spectrum  is  Ignored  and  the  next  spectrum  is  analyzed.  If  the  diagonal  element  is  incremented 
fay  one.  we  also  Increment  by  one  all  elements  of  the  matrix  for  the  row  and  column 
containing  the  diagonal  element.  The  elements  of  the  row  are  incremented  because  If  we  see 
an  event,  for  a  given  pass,  that  satisfied  the  threshold  requirements  In  one  latitude  zone,  then 
starting  at  that  latitude  and  Integrating  spatially  either  equatorward  or  poleward  the 
probability  of  seeing  an  event  for  that  pass  must  be  unity.  Similarly,  the  column  is 
Incremented  because  starting  at  any  other  latitude  for  that  pass  and  Integrating  either 
poleward  or  equatorward.  the  probability  must  be  unity  If  the  Integral  extends  to  the  latitude 
where  the  event  satisfying  the  thresholds  was  observed.  The  spectra  taken  within  one  auroral 
pass  will  contribute  to  a  number  of  dlfierent  arrays  and  one  spectrum  can  contribute  to  more 
than  one  array.  For  each  pass,  a  counter  for  each  MLT  sector  in  each  array  is  Incremented  if. 
for  that  pass,  that  satellite  passed  through  that  MLT  sector.  The  desired  probabilities  are 
calculated  by  dividing  the  values  In  the  matrix  by  the  number  of  ptasses  occurring  In  that  MLT 
sector  . 

A  simplified  case  is  used  to  illustrate  the  procedure.  A  simplified  latitudinal  profile  of  the 
Integral  energy  flux  is  shown  in  Figure  2.  For  the  purpose  of  this  demonstration  we  assume 
that  all  the  electron  spectra  in  the  pass  have  an  average  energy  of  greater  than  1  keV.  that  the 
satellite  is  In  the  2100-2400  MLT  sector  and  that  Kp  is  In  the  range  fi'om  2-  to  3+.  For  this 
pass  we  will  determine  the  inputs  to  the  matrix  for  the  Integral  energy  flux  greater  than 
3.16x10®  keV/(cm^  sec  ster).  In  the  figure,  vertical  lines  are  drawn  at  the  2®  intervals 
corresponding  to  each  of  the  diagonal  elements.  A  horizontal  line  Is  drawn  representing  the 
threshold  for  the  matrix.  From  the  figure  we  determine  that  for  each  2°  element  In  latitude 
only  between  62  and  70®  CGL  Is  there  at  least  1  second  in  each  element  when  the  Integral 
energy  flux  was  above  the  threshold,  'fhe  five  diagonal  elements  for  this  latitude  range  are 
Incremented  by  one  and  the  rows  and  columns  are  filled  in  as  described  m  the  previous 
example.  This  is  shown  In  Table  2. 

For  such  a  simplified  case  we  can  test  whether  the  matrix  produces  the  correct  Intv^ral 
probabilities.  If  we  ask.  "what  Is  the  probability  of  seeing  an  event  above  threshold  at  50® 
CGL?".  It  corresponds  to  the  matrix  element  in  the  extreme  lower-left  comer  in  Table  2.  This 
Is  zero.  In  agreement  with  Figure  2.  To  determine  the  Integral  probability  as  we  move  to 
higher  latitude  starting  at  50®  CGL  we  must  consider  the  lowest  row.  This  Integral  probability 
Is  zero  until  the  sixth  column  and  unity  thereafier.  This  Is  again  consistent  with  the  figure. 
From  50®  CGL  poleward,  no  event  above  the  threshold  is  seen  at  latitudes  less  than  60®.  This 
corresponds  to  the  first  five  columns  ^shlch  are  indeed  zero.  Since  the  threshold  Is  exceeded 
between  60  and  62®  CGL,  the  integral  probability  from  50®  to  any  latitude  above  60®  will  be 
unity  as  is  the  case  In  the  matrix. 

This  same  procedure  was  repeated  for  all  auroral  passes  of  the  F9  DMSP  satellite  occurring 
In  1990.  This  represented  10,141  auroral  passes  and  a  total  of  over  12  million  individual  1- 
second  sp>ectra.  This  data  set  was  of  sufficient  size  that  a  good  statistical  determination  was 


50 
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80 
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Figure  2.  Energy  Flux  (JEJtot^  Spectrum  for  a  1/4  Satellite  Orbit  (or  1/2 
Hemispheric  Pass).  JEijQt  is  given  in  KeV/(cm2  sec  ster)  as  a  function  of 
corrected  geomagnetic  latitude  (CGL). 
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Table  2.  Probability  Matrix  (i.j)  for  One  Pass  Through  the  Auroral  Oval 


20  1111111 
19  1111111 

18  1111111 
17  1111111 

16  1111111 
15  1111111 

14  1111111 

13  1111111 

12  1111111 
U  1111111 
IQ  1111111 
9  1111111 

8  1111111 
7  1111111 

6  1111111 
5  0  0  0  0  0  1  1 
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obtained  of  both  the  probability  of  occurrence  above  the  threahold  at  Individual  latitude  zones 
with  the  MLT  coverage  of  the  satellite  and  also  for  the  probability  Integrated  over  latitude. 


3.  RESULTS 


Figure  3  shows  a  typical  example  of  the  output  of  the  stucty.  This  figure  Is  generated  using 
the  software  package  EPDF9.  A  step-fay-step  guide  to  the  operation  of  this  software  Is  provided 
In  the  Appendix.  The  case  shown  Is  for  an  average  energy  greater  than  7  keV.  an  Integral 
energy  flux  greater  than  50  eig/(cm^  sec)  (or  50  kR).  far  the  2100-2400  MLT  sector,  and  for 
In  the  range  from  2-  to  3+.  The  rectangular  section  is  a  color  spectrogram  of  the  results.  The  y 
axis  gives  the  delta  latitude  equatorward  or  poleward  of  a  specific  latitude  over  which  the 
probability  Is  integrated.  The  x  axis  gives  tl.c  latitudes  fix>m  50  to  90^  CGL  wlthm  the  specified 
MLT  sector.  Along  the  zero  delta  latitude  line,  the  color  scale  gives  the  probability  of 
observing  an  event  above  the  threshold  In  each  2°  latitude  zone.  For  each  CG  latitude,  the 
color  scaling  above  and  below  the  zero  delta  latitude  line  gives  the  probability  of  seeing  an 
event  above  those  thresholds  Integrated  over  a  ddta  latitude  poleward  (positive  delta  latitude) 
or  equatorward  (negative  delta  latitude)  starting  at  that  CG  latitude.  The  bar  graph  at  the 
bottom  gives  the  number  of  passes  in  which  an  event  was  seen  in  each  of  the  20  CGL  zones. 

For  this  example,  one  sees  that  the  maximum  probability  of  seeing  an  event  within  a 
sln^e  CGL  zone  occurs  between  66  and  68^  CGL  (approximately  1 .4  percent).  Integrating 
poleward  of  that  latitude  over  a  delta  latitude  of  8**  Increases  that  probability 
approximately  a  factor  of  3.  Similarly,  inte^tlng  equatorward  of  that  latitude  over  a  delta  of 
6°  approximately  doubles  the  probability.  The  minimum  spatial  mtegratlon  required  to 
obtain  the  highest  probability  of  approximately  5.6  percent  Is  approxlmatety  30"  CGL  starting 
at  58"  CGL  and  Integrating  poleward  or  25"  starting  at  75"  and  integrating  equatorward. 

We  next  consider  how  the  maximum  and  Integral  probabilities  vary  as  a  function  of  the 
variation  of  thresholds  and  spatial  and  activity  separations.  We  first  examine  the  variation 
with  changes  In  the  mtegral  energy  flux  threshold.  For  this  case,  we  take  an  average  energy 
threshold  of  7  keV  for  the  2100-2400  MLT  sector  and  for  the  Kp  range  finm  2-  to  3+.  The 

a 

results  are  summarized  In  Table  3  and  are  shown  as  gray  scale  plots  In  Figures  4  through  9.  In 
the  table,  the  left  to  the  right  colunms  are  the  energy  flux  thresholds,  the  maximum 
probability  seen  In  a  single  2"  latitude  zone,  the  corrected  geomagnetic  latitude  of  the 
maximum,  vdiether  a  6"  latitude  mtegratlon  from  the  maximum  probability  gives  the  largest 
value  for  mtegratlng  poleward  (N)  or  equatorward  (S).  the  probability  this  mt^ratlon  gives, 
and  the  largest  mt^ral  probability  possible. 

For  this  case,  at  lower  mtegral  energy  flux  thresholds,  the  maximum  probability  Is  fotmd 
at  a  latitude  of  65"  CGL  and  moves  to  67"  for  the  50  and  160  kR  thresholds.  This  probability 
decreases  with  mcreaslng  mtegral  energy  flux  thresholds,  dropping  from  9  percent  at  the 
lowest  threshold  to  0.3  percent  at  the  highest.  Startmg  at  the  maadmum  probability  pomt  and 
mtegratlng  northward  over  a  delta  latitude  of  6".  these  probabilities  can  be  significantty 
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Figure  3.  A  Typical  Plot  of  the  Integral  Probability  of  Auroral  Electron 
Flux  Barents.  The  parameter  ranges  are  provided  in  the  legend. 
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Events 


Table  3.  Maximum  Probability  as  a  Function  of 
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Maximum 
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12.6 

16.2 
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160 

0.3 

67 

N 

1.2 

2.7 

CGL 


Per-Pass 

Percentage 

Probability 

Eave  (KeV) 
>=  7 

JEtot  (kR) 
>=  0.5 


MLT 

21  -  24 


1.8 

.0 


Figure  4.  Integral  Probability  of  Auroral  electron  Flux  Events 
Greater  or  Equal  to  0.5  kR  (or  1.0  erg/sec  ster).  Other  parameter 
ranges  are  provided  in  the  legend. 
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Figure  6.  Same  as  Figure  4  Except  that  the  Energy  Flux  JE^tot  ^ 
5.0  kR  (or  10.0  erg/sec  ster). 
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F^re  8.  Same  as  Figure  4  Except  that  the  Energy  Flux  JE^ot  ^ 
50.0  kR  (or  100.0  erg/sec  ster). 
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Figure  9.  Same  as  Figure  4  Ebccept  that  the  Energy  Flux  JE^ot  ^ 
160.0  kR  (or  320.0  eig/sec  ster). 
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Increased.  For  the  three  lowest  thresholds  the  Increase  is  approximately  50  percent.  For  the 
higher  thresholds  the  Increase  can  be  as  much  as  a  factor  of  4.  The  maximum  spatially 
Integrated  probabilities  are  even  larger,  varying  from  16.2  percent  at  the  lowest  threshold 
(approximately  a  factor  of  2  Increase)  to  2.7  percent  at  the  highest  threshold  (approximately  a 
factor  of  9  Increase). 

We  next  consider  the  case  where  we  vary  the  average  energy  threshold.  For  this  case  we 
choose  an  Integral  energy  flux  threshold  of  5  kR  for  the  Kp  range  of  2-  to  3+  and  for  the  MLT 
sectors  2100-2400  MLT  and  0600-0900  MLT.  The  results  are  summarized  In  Table  4.  For  both 
local  time  sectors  the  maximum  probability  decreases  with  Increasing  average  energy,  falling 
from  12  to  0.2  percent  for  the  0600-0900  MLT  sector  and  from  20  to  0.2  percent  for  the  2100- 
2400  MLT  sector.  For  the  0600-0900  MLT  sector  the  location  of  the  maximum  moves  to  a 
lower  latitude  with  Increasing  average  energy  while  for  the  2100-2400  MLT  sector,  the 
maximum  Is  located  at  65®  CGL  for  all  but  the  largest  average  energy.  As  for  the  previous  case, 
the  probabilities  are  significantly  Increased  Integrating  over  latitude.  Integrating  over  6° 
from  the  maximum  Increases  the  probability  by  50  to  60  percent  at  the  lower  three  thresholds 
for  both  MLT  sectors  and  by  more  than  100  percent  at  the  highest  thresholds. 

Table  4.  Maximum  Probability  as  a  Function  of  Ejive 
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Threshold 

Maximum 
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Probability 
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Direction 
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1.0 
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65 

N 

32 

36 

3 

14 

65 

N 

24.5 

28 

7 

6 

65 

N 

12.8 

14.4 

10 

3.2 

65 

N 

4.8 

13. 

13 

1.05 

65 

N 

2.1 

3.5 

17 

0.2 

67 

N 

0.8 

1.6 

Lastly,  we  consider  the  case  where  the  thresholds  and  MLT  sector  are  held  constant  and  the 
Kp  range  is  varied.  For  this  case  the  integral  energy  flux  threshold  is  5  kR,  the  average  energy 
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threshold  Is  7  keV  and  we  looked  at  the  2100-2400  MLT  and  0600-0900  MLT  sectors.  These 
results  are  summarized  in  Table  5.  One  notes  that  for  these  thresholds  the  maximum 
probabilities  for  both  local  times  occur  for  Kp  in  the  range  from  2-  to  3+.  The  values  of  the 
maximum  probability  for  the  two  cases  are  also  approximately  the  same  for  the  MLT  sectors. 
The  location  of  the  maximum  probability  moves  equatorward  in  latitude  with  increasing 
activity  In  both  MLTs;  from  67  to  63°  for  the  0600-0900  MLT  sector  and  from  67  to  61°  for  the 
2100-2400  MLT  sector.  This  behavior  Is  consistent  with  the  known  expansion  of  the  oval  to 
lower  latitudes  with  increasing  activity.  The  Increase  In  the  probability  with  the  Integration 
In  latitude  Is  much  the  same  as  for  the  previous  cases  shown. 


Table  5.  Maximum  Probability  as  a  Function  of  Kp 


Kp  range 

Maximum 

Probability 

CGL 

North/ South 
Direction 

±  6°  ACGL 
Probability 

Maximum 

Integral 

Probability 

0600-0900  MLT 

0.7 

67 

S 

1.4 

3.15 

2-  -»3+ 

9 

67 

N 

12.6 

16.2 

4-  -^5+ 

4.5 

65 

N 

7.2 

8.1 

6 — >  9q 

0.6 

63 

N 

1 

1 

2100-2400  MLT 

Oo  ->  1+ 

7 

67 

N/S 

2.10 

3.15 

2- -^3+ 

8 

65 

N 

12.8 

14.4 

4-  -»5+ 

4.5 

63 

N/S 

6.3 

8.1 

6-  -*  9q 

0.75 

61 

N/S 

1.2 

1.5 

4.  CONCLUSION 


A  statistical  study  has  been  completed  using  the  SSJ/4  data  from  the  F9  DMSP  satellite. 
The  results  of  the  study  provide  the  capability  to  determine,  for  locations  about  the  auroral 
oval,  the  probability  of  observing  auroral  electron  precipitation  above  specified  levels  In  the 
Integral  energy  flux  and  average  energy  for  any  line  of  sight  through  the  oval.  This  provides  a 
new  tool  for  satellite  flight  mission  planners  to  determine  the  pointing  of  onboard  ^sterns  to 
maximize  the  probability  of  observing  specific  levels  and  types  of  auroral  phenomena.  More 
Importantly,  the  results  provide  the  capability  to  determine  for  operational  Air  Force  systems 
the  probability  that  auroral  precipitation  at  specific  level  Is  producing  backgroimd  to 
surveillance  systems  operating  through  the  oval  or  Ionospheric  perturbations  for  C^I 
operating  In  and  through  the  auroral  oval. 


17 


The  only  deflctencies  In  the  present  study  arises  from  the  use  of  only  one  year  of  DMSP  F9 
SSJ/4  data.  The  use  of  only  the  SSJ/4  data  fixjm  the  F9  satellite  limits  the  magnetic  local 
time  coverage  of  the  study  and  the  use  of  only  one  year  limits  the  statistical  accuracy  of  the 
results.  These  {problems  can  be  eliminated  through  the  use  of  data  from  the  SSJ/4  detectors  on 
the  F7,  F8.  FIO  and  Fll  satellites.  In  combination,  these  instruments  have  recorded  many 
years  of  data  and  provide  nearly  complete  MLT  coverage  of  the  oval.  Such  an  extension  of  the 
present  study  remains  an  option. 
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Appendix 

Description  of  the  PC  Software  EPDF9 


The  software  that  generates  the  integral  probability  plots  is  interactive  and  is  menu-driven 
using  the  binned  data  base  containing  the  total  number  of  auroral  event  occurrences.  This 
software  was  tested  on  a  IBM  PC.  The  graphic  development  was  based  on  a  commercial  graphics 
package  (Graphic  by  Scientific  Endeavors  Inc.).  Installing  the  software  in  other  PC  machmes  is  a 
simple  procedure.  The  EPDF9  software  uses  a  menu  with  several  options  given  in  Table  A1 .  The 
first  four  options  set  the  auroral  peuameter  ranges.  These  parameters  are  the  eneigy  flux 
threshold,  the  average  energy  threshold,  the  magnetic  local  time  sector,  and  the  Kp  range.  Option 
5  initiates  the  plot  once  the  parameters  have  been  selected.  Option  6  is  to  save  the  plot  (*.TKF)  and 
listing  (*.LST)  files.  Option  7  makes  it  possible  to  display  the  results  on  the  screen  or  generate  a 
laser  copy.  Option  8  exits  the  program.  Option  9  indicates  that  the  user  has  finished  selecting  the 
p2U’ameters. 

Teible  A2  lists  the  possible  inputs  for  each  option.  The  value  to  be  given  as  the  input  is  the 
number  in  parentheses  for  each  of  the  first  four  options.  The  user  is  eisked  to  change  the 
peirameter  values  of  the  options  that  were  selected  in  the  previous  menu.  The  program  interprets 
all  options  in  the  order  selected  or  until  a  plot  is  requested.  The  plots  are  then  generated  on 
request.  The  program  again  displays  the  menu  in  Table  A 1  if  no  plot  is  requested  or  after  each 
plot  is  completed.  This  procedure  is  repeated  until  the  user  selects  Option  8,  which  ends  the 
session. 

The  results  are  displayed  in  terms  of  int^ral  probability  versus  latitude  zones  for  a  fixed 
local  time  sector,  average  eneigy  and  energy  flux  thresholds.  The  probabilities  are  referred  to  as 
the  per-pass  probability  since  they  are  derived  from  measurements  collected  during  individual 
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Individual  hemispheric  passes.  The  number  of  events  observed  In  each  latitude  zone  Is  also 
shown  In  a  lilstogram  In  the  bottom  of  each  plot.  This  histogram  should  Indicate  how  reliable 
the  computed  probabilities  are.  The  gray  scale  plots  were  generated  for  this  document  on  the 
HP  LaserJet  m. 


Table  Al.  EPDF9  Menu 


1)  JEtot 

2)  Eave 

3)  Local  Time 

4)  IQ)  Range 

5)  Start  Plot 

6)  Save  Files 

7)  Screen/Printer 

8)  Quit  Session 

9)  Done  with  Selection 

Table  A2.  EPDF9  Option  Description 


Option 

Description 

1 

Energy  Flux  Thresholds  In  kR 
(1)  0.5(2)  1.6(3)  5.0 
(4)  16.0  (5)  50.0  (6)  160.0 

2 

Average  Energy  Thresholds  In  keV 
(1)  1.0  (2)  3.0  (3)  7.0 
(4)  10.0  (5)  13.0  (6)  17.0 

3 

local  time  3-hour  sector 

(1)  0  -  3  (2)  3  -  6  (3)  6  -  9  (4)  9  -  12 

(5)  12  -  15  (6)  15  -  18  (7)  18  -  21  (8)  21  -  24 

4 

IQ)  Index  range 
(1)  Oq  .  1+  (2)  2-  .  3+ 

(3)  4-  .  5+  (4)  6-  .  9o 

5 

Generate  plot  and  listing 

6 

Save  plot  (.TKF)  and  listing  (.LST)  files 

7 

Screen  (0)  or  Laser  Printer  (1)  Option 

8 

Terminate  Session 

9 

Quit  Menu 

We  next  Illustrate  a  typical  sample  run.  To  be  able  to  run  EPDF9  on  a  PC  system,  the 
following  files  must  be  available  In  the  default  directory: 
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EPDF9.EXE  main  executable  program 

MSXDBA.DAT  binned  DMSP  data  file 

MSXPSD.DAT  total  number  of  passes  In  each  spatial  grid 

*.FNT  graphics  font  files 

*.01^  graphics  screen  files 

*.PRN  graphics  printer  files 

SAVER.BAT  batch  file  to  save  plot  and  listing  files 

The  following  steps  Illustrate  a  sample  run  of  the  EPDF9  software. 

Step  1:  ‘lype  EPDF9  and  a  menu  Is  displayed  as  In  Table  Al.  We  choose  the 

following  options  for  this  session: 

132459 

These  options  will  be  executed  in  the  order  given.  The  first  four 
options  request  changes  of  the  energy  flux  JE^ot  threshold,  the  MLT 
sector,  the  average  energy  E^^  threshold,  and  the  Kp  range,  In  this 
order.  Option  5  Is  then  requested  for  generating  a  screen  plot 
(^ilch  Is  the  default  at  the  beginning  of  each  session).  Option  9  Is 
terminate  the  Input  selections. 

Step  2:  Option  1:  To  change  JE^^^  threshold,  use  one  of  the  values  In  Table 

A2.  We  choose  bln  #  4;  te.,  all  energy  fluxes  greater  or  equal  to 
16  kR.  Note  tliat  to  keep  the  current  JE^^j  range,  enter  any 

n^atlve  integer  value  can  be  entered. 

Option  3:  To  change  MLT  sector,  use  one  of  the  values  In  Table  A2. 
We  choose  bln  #  8;  te.,  all  measurements  that  were  recorded  In  the 
2100  to  2400  MLT  sector.  Again  to  keep  the  current  MLT  sector, 
enter  any  negative  integer  value. 

Option  2:  To  change  the  average  energy  E^^  threshold,  use  one  of 

the  values  In  Table  A2.  Use  the  defiiult  value  (bln  #  1)  by  entering 
any  negative  Integer.  This  means  that  all  measurements  with  E^^ 

greater  or  equal  to  1  keV  are  considered.  Note  that  option  2  comes 
after  option  3  due  to  the  selection  order  In  step  1. 


Option  4:  To  change  the  gecanagnetic  activity  Kp  range,  use  one  of  the 
ranges  tai  Table  A2.  In  this  case,  two  values  must  be  entered.  We 
choose  the  lower  bound  and  the  upper  bound  of  Kp  to  be  In  bln  #  2;  Le.. 
all  measurements  that  were  recorded  when  Kp  was  between  2-  and  3+. 

Option  5:  This  option  generates  a  plot.  Since  option  7  was  not 
selected,  a  plot  on  the  screen  will  be  generated,  which  is  a  default  at 
the  beginning  of  each  session. 

At  this  point,  the  data  is  read  In  and  the  appropriate  steps  are  taken 
to  generate  the  percent  probabilities  for  the  specified  electron  event. 

Step  3;  After  examination  of  the  plot,  pressing  the  RETURN/ENTER  key 

displays  the  menu  in  Table  A1  as  in  Step  1.  We  may  select  Option  6  to 
save  the  listing  and  plot  files  or  generate  an  HP  laser  copy  of  this  plot. 
To  generate  a  hardcopy,  we  enter  the  7.  5.  9  choices,  in  this  order. 

Step  4:  Option  7:  To  obtain  a  hardcopy  of  this  plot,  we  choose  the  printer 

option  by  entering  1. 

Option  5:  Generates  the  plot  and  sends  the  results  to  the  HP  laser 
printer.  This  plot  is  shown  in  Figure  Al.  The  Integral  probability 
versus  latitude  is  given  on  the  top  panel  of  Figure  Al  for  this  electron 
event.  The  cut-off  of  each  probability  range  is  given  by  a  shaded 
pattern  on  the  right-hand  side  of  the  figure.  The  horizontal  axis 
shows  the  corrected  geomagnetic  latitude  (CGL)  in  degrees.  The 
vertical  axis  gives  the  offset  in  degrees  from  CGL;  in  other  words,  the 
vertical  axis  is  the  delta  latitude  (ACGL).  When  ACGL  is  not  zero,  the 
probability  vedue  is  the  cumulative  probability  from  ACGL*0  to  this 
latitude  offset.  Positive  ACGL  means  that  we  are  looking  poleward 
from  CGL.  Negative  ACGL  means  that  we  eae  looking  equatorward 
from  CGL.  The  bottom  panel  shows  the  number  of  occurrences  of  the 
electron  events  over  all  latitude  zones. 

Step  5;  After  the  plot  is  generated,  pressing  the  RETURN/ENTER  key  again 

displays  the  menu  in  Table  Al.  To  terminate  the  session,  we  type  in  8. 
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Figure  Al.  Sample  Plot  of  the  Integral  Probability  of  the  Auroral 
Precipitating  Electron  Events  with  the  Parameter  Ranges  Provided  in 
the  Legend. 


